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Recent experiments on microscopic graphite mesas demonstrate reproducible high-speed microscale 
superlubricity, even under ambient conditions. Here, we explore the same phenomenon on the nanoscale, by 
studying a graphene flake sliding on a graphite substrate, using molecular dynamics. We show that 
superlubricity is punctuated by high-friction transients as the flake rotates through successive crystallographic 
alignments with the substrate. Further, we introduce two novel routes to suppress fi-ictional scattering and 
acliieve persistent superlubricity. We use graphitic nanoribbons to eliminate frictional scattering by 
constraining the flake rotation, an approach we call frictional waveguides. We can also effectively suppress 
frictional scattering by biaxial stretching of the graphitic substrate. These new routes to persistent 
superlubricity at the nanoscale may guide the design of ultra-low dissipation nanomechanical devices. 

The nonlinear nature of friction at the nanoscale has sparked intense experimental and theoretical interest^l 
A focus of much of this research is on achieving superlubricity^"^, a term originally introduced to describe 
ultra-low friction between rotationally misaligned - and hence incommensurate - crystalline lattices^"^. 
Superlubricity has been observed by Friction Force Microscopy, which typically operates at ^1 mm/s. However, 
for freely-moving nanoparticles at low velocities, superlubricity is suppressed because of spontaneous rotation 
into energetically favourable commensurate states, or "lock-in"^°'^\ while the commensurate states are always the 
stable orbits for the sliding of the graphene flake^^'^^. Recent experimental results, however, show reproducible 
high-speed superlubricity in microscale graphite mesas (10 X 10 ^m^) under ambient conditions, where the 
highest speed is 25 m/s^^. Besides, no lock-in is observed for a wide-range of rotational misalignments, due to 
much larger rotational inertia on the microscale^^. In this study, we wish to explore in detail how the lock-in 
phenomenon disrupts superlubricity on the nanoscale, and consequently how we may design systems that 
suppress this effect, so that persistent superlubricity can be exploited for building efficient nanomechanical devices. 

In recent simulations by Guerra et al/^, a regime of "ballistic nanofriction" was observed for nanoclusters of 
gold moving at high velocity on graphite. A profound insight of this study is that in contrast to the much- studied 
low- velocity regime, at high velocities superlubricity is the norm, punctuated by transitory periods of lock- in 
where friction is high. Lebdeva et al.^^ have observed similar behavior in simulations of the thermally induced 
motion of nanoscale graphene islands on graphite. 

Results 

Here we introduce a simple mechanism to induce ballistic nanofriction in graphene nanoparticles, based on high- 
speed motion observed experimentally by shearing microscopic mesas of graphite, which occurs in order to 
minimize the free surface energy at the sheared interface^^'^^'^^. From classical mechanics, it is straightforward to 
deduce an expression for the peak speed produced by this effect for shearing of a square flake of graphite from a 
fixed mesa: v^ax = [(2/iV)(s//)(£f/mc)] where AT is the number of graphitic layers in the flake, £f is the exfoliation 
energy per carbon atom and rric is the mass of a carbon atom, s is shear displacement of the graphite flake and / is 
the length of the flake. Recent experiments have shown the microscopic graphite flakes of thickness about 100 nm 
can reach 25 m/s, in quantitative agreement with the above formula^^. For graphene, by extrapolating toN= 1, we 
deduce retraction speeds >400 m/s should be attainable. 
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Figure 1 | Retraction motion of the graphene flake. Four snapshots (a) of 
the graphene flake of initial commensurate state during retraction where 
the direction and magnitude of the shding speed are indicated by the 
arrows. The retraction speed and retraction force in x direction for initial 
commensurate and incommensurate state are recorded in (b). The four 
numbers are corresponding to the four snapshots in (a), (c) The rotation 
angle and angular velocity of the flake. 

To investigate the sliding behavior of a two-dimensional nanoscale 
graphene flake at such speeds, we employ molecular dynamics simu- 
lations using the Gromacs code^^ with Drieding force field to describe 
both the intra-layer and inter-layer interaction between carbon 
atoms of graphene^^. The intra-layer interaction is consisted by bond 
stretching, bond-angle bending and dihedral angle torsion, whfle the 
inter-layer interaction is described by Lennard- Jones van der Waals 
potential. The parameters in the potential function are validated by 
correctly characterizing the structural and mechanical properties of 
different carbon nanostructures^^"^^. A graphite substrate is modeled 
by two layers of graphene in A-B stacking with the size of 30 nm X 
20 nm where the bottom layer is fixed during the simulation. The 
moving graphene flake is modeled by one layer graphene with the 
size of 10 nm X 10 nm. Here two types of motion are studied in our 
simulations: the retraction motion and high speed sliding of the 
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Figure 2 | Flake dynamics after launch at 400 m/s. Four snapshots (a) of 
the instantaneous orientation where the arrows indicate the direction and 
magnitude of the sliding speed and the frictional scattering at four points 
indicated on the flake trajectory (b). The speed components of the flake 
and Vy (c), show stepwise changes at each scattering. The flake rotation 6 
(d) shows frictional scattering occurring primarily at high-symmetry 
angles (multiples of 60°). The force components acting on the flake, and 
Fy (e) show large oscillatory variations at each scattering. The kinetic 
energy component decays within about 1 ns to less than 10% of its 
initial value (f), and becomes comparable to Ey and E^. 

graphene flake, where the time step in ah of our simulations is Ifs. 
For the retraction motion the top graphene flake is sheared away 
from the substrate with a half of its length s = 1/2, as iflustrated in 
Fig. la and free boundary conditions are applied to the substrate. For 
the simulation of the high speed sliding the graphene flake is placed 
on the center of the substrate and periodic boundary conditions are 
applied to the substrate. The whole system is first energy minimized 
by conjugate gradient method, and then the sliding speed is applied 
to the graphene flake. Here, the micro-canonical assemble is used in 
ah of the simulations in coordination with the previous works of 
studying the frictional behavior of graphene^^'^^. 

The suspended graphene flake is released to simulate the retrac- 
tion motion and the retraction speed, force, rotation angle and angu- 
lar velocity of the graphene flake are recorded, as shown in Fig. Ib-c. 
Here, the retraction speed is the speed of the center of mass (COM) of 
the sliding flake; the rotation angle is defined as the angle rotated to 
its COM. The retraction force is calculated as the force between the 
substrate and the flake during retraction through the Lennard- Jones 
van der Waals potential defined before. Two types of initial condi- 
tions are considered: the first is the graphene flake commensurate 
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Figure 3 | Close-up of frictional scattering. Variation of the speed component (a) and the van der Waals bonding energy £vdw (b) as flake displacement 
along the x-direction, before, during and after a scattering event. Numbers indicate the snapshots of the flake (c). Vertical dotted lines in (a) and (b) 
indicate the angular width AO^ of the coherent region where a^> 21. The instantaneous corrugation potential, E^y felt by the flake, for relative 
displacement in the x-direction about the flake's actual position (indicated by a circle) is shown before (d), during (e,f) and after (g) frictional scattering. 



with the substrate; the second is the graphene flake incommensurate 
with the substrate where the orientation of the moving graphene 
flake is 90° misaligned to the substrate. Although the suspended 
graphene flake is very easy to fold and wrinkle at finite temperature, 
the length of the suspended graphene flake in our simulation is 5 nm 
which is hard to fold. Besides, the initial temperature of the whole 
system is 0 K and after releasing the graphene flake quickly retracts 
back, so the flat configuration is well kept. The maximum retracting 
speeds for the two cases are about 600 m/s which are consistent with 
the simple formula prediction. However, for the commensurate case 
the retraction is accompanied by the rotation at very beginning of the 
retracting which is corresponding to a large fluctuation of the retrac- 
tion force, see the retraction force in Fig. lb and angular velocity in 
Fig. Ic. After that the angular velocity is almost constant and the 



fluctuation of the retraction force quickly decreases. As the rotation 
the final retraction speed for the commensurate state is smaller than 
the incommensurate state. 

Based on our results the commensurate state seems unstable for 
high sliding speed of the graphene flake and the graphene flake will 
automatically rotate to avoid the commensurate state. This is differ- 
ent from the previous works at low sliding speed where the commen- 
surate state is the most stable state^^'^^. In what foflowing we will 
simulate the high-speed sliding behavior of the graphene flake and 
find ways to reduce the sliding friction in such speed. 

Results presented in Fig. 2 are for a 10 X 10 nm^ fiake launched at 
^max — 400 m/s along the armchair direction of the substrate by 
adding the same velocity 400 m/s in x direction to every atom in 
the sliding flake. The sliding distance corresponding to 1 time step 
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(Ifs) at Vmax is Icss than 0.0004 nm. The initial state of the flake and 
substrate is an energy minimum state with A-B stacking ahgnment, 
and the initial temperature of the whole system is 0 K. The sliding 
speed, rotation angle, angular velocity and force are defined same as 
the ones in the retraction simulations. The sliding behavior in 
Figure 2b shows the trajectory of the flake, with sporadic changes 
in the direction of motion indicated by vertical dashed lines. For 
reasons that will become apparent shortly, we refer to these events 
as frictional scattering. Snapshots in Fig. 2a show each scattering 
event occurs when the flake is briefly crystallographically aligned 
with the substrate. At lower speed, this alignment would result in 
lock-in. But at high speed, the flake has sufficient kinetic energy to 
escape lock-in. 

Figure 2c shows the component of the sliding speed along the x- 
direction, v^, dropping in discrete steps at each scattering, with sim- 
ilar stepwise variations for Vy. The evolution of flake rotation angle, 
plotted in Fig. 2d, indicates that the flake angular velocity, o), also 
changes abruptly at each scattering which is similar to the rotation in 
the retraction in Fig. Ic. Fig. 2e shows that scattering events subject 
the flake to large oscillatory forces in both x and y directions, pro- 
ducing abrupt changes of direction and velocity. In Fig. 2f, the linear 
kinetic energy components and Ey in the x- and )/- directions and 
the rotational kinetic energy of the flake are displayed. For better 
comparison we enlarge 10 times of Ey and E(^ in Fig. 2f Since Ey and 
E^ remain small throughout, we conclude that the initial kinetic 
energy of the flake is primarily dissipated into vibrational energy, 
with enhanced dissipation at each scattering. After ~ 1 ns, the sliding 
speed is smaller than 100 m/s, with E^ and Ey reaching similar mag- 
nitudes, and the flake moving randomly. Since the micro -canonical 
assemble is applied in the simulation, the whole system will be heated 
up as the dissipation of the translational speed of the sliding flake. 
However the substrate in our simulation (30 X 20 nm) is much 
larger than the sliding flake (10 X 10 nm), we estimate the temper- 
ature raising after dissipating all of the translational kinetic of the 
sliding flake is smaller 10 K. 

To study the effect more closely, we zoom in on a region of fric- 
tional scattering in Fig. 3. The flake speed component, v^, and van der 
Waals bonding energy between flake and substrate, £vdw which is the 
van der Waals energy between the flake and substrate during the 
sliding, are displayed in Figs. 3a-b, respectively, and show rapid 
oscillations near alignment. Fig. 3c shows eight snapshots of the 
moire pattern generated by the overlap of flake and substrate lattices, 
at different points indicated in Figs. 3a-b. When the flake is rota- 
tionally misaligned with the substrate, the moire pattern has a small 
unit cell size, a^. This means the spatial phase of flake atoms relative 
to the substrate potential oscillates rapidly across the flake. The result 
is a very small net corrugation potential for the flake as a whole: this is 
the origin of superlubricity for incommensurate lattices. 

As the flake rotates into alignment, expands (snapshots 1-4). 
Near alignment, once ^ /, where / is the linear size of the flake, 
all flake atoms are nearly in phase with the substrate (snapshots 5, 
6) and the net corrugation potential experienced by the flake 
increases abruptly, resulting in scattering. The corrugation poten- 
tial drops again once the flake rotates sufficiently far out of align- 
ment (snapshots 7, 8). Here the corrugation potential E^ is defined 
as the van der Waals energy between the flake and substrate by 
displacing the flake relative to flake's actual position along its 
instantaneous sliding speed, and recalculating £vdw keeping all 
other parameters flxed. The variation of the corrugation potential, 
£cj before, during and after alignment is shown in Figs. 3d-g. Near 
alignment (snapshots 5, 6) E^ is —lOOx larger than during super- 
lubricity, and the flake is subject to large potential gradients as 
indicated in Figs. 3e, f. 

In the reference frame of the moving flake, the fundamental 
Fourier component of the substrate potential along the x-direction 
appears as a travelling wave of wavelength X = d, where d is the in- 



plane lattice constant of graphite. When rows of atoms in the flake 
align with crests of this wave, there is a maximum oscillatory force 
acting on the flake. This corresponds to the condition for Bragg 
scattering IdsmO = nk, at normal incidence with n = 2. (For odd 
n, the force on neighbour rows cancels out.) 

Since frictional scattering depends on the crystallinity of the flake, 
it will not be observed for diffusion of single atoms or small molecules 
at similar velocities'^. Also, since frictional scattering depends on 
rotation in a plane, it is an intrinsically two-dimensional effect, phys- 
ically distinct from velocity- dependent resonances observed in one- 
dimensional models^^ and nanotube simulations'\ 

We can derive the angular width JOc of the coherent region where 
a^i ^ / using the standard moire fringe formula'^ = d/(2\sin(J0/ 
2) I ), where AO is the rotation angle relative to alignment. Substituting 
for = 2/ as a lower bound for coherence across the flake gives AO 
= d/2l in the small angle limit. Since AO is symmetric around the 
aligned position, AOc = d/l. This is identical to the Scherrer equation 
for a crystallite of dimension / and radiation of wavelength X = d, for 
the case of normal incidence, emphasizing the mathematical analogy 
with Bragg scattering, which has been employed in the previous 
works to describe the interlocking between two contact surfaces'^'^^. 
Substituting d = 0.246 nm, the lattice constant of graphite and / = 
10 nm gives AOc= 1.4°. This angular width is plotted in Fig. 3b, and 
agrees with the range where £vdw oscillates rapidly. 

How robust an effect is frictional scattering? We have studied it as 
a function of several parameters, such as size of the graphene flake, 
temperature of the system and initial orientation of the graphene 
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Figure 4 | The influence of flake size. Comparison of decay of the speed 
component versus time after launch for different flake sizes (a) shows 
similar behavior for all sizes, (b) The duration of intervals of superlubricity, 
Ts, depends inversely on the average angular velocity of the flakes during 
the intervals. 
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Figure 5 | The sliding behavior for different initial conditions, (a) The 

sliding behavior of 5 X 5 nm^and 10 X 10 nm^ flakes launched initially at 
0 K and 300 K. (b) The shding behavior of 10 X 10 nm" flake for different 
orientation (the misalignment between the flake and the substrate is from 
0° to 30°). 

flake. For example, we find that for flakes with sizes from 4X4 nm^ 
to 10 X 10 nm^, the overall behavior is simflar; the forward sliding for 
all flakes being dissipated within 2 ns, see Fig. 4a. The sliding speed in 
X direction sometimes increases when it is smaller than 100 m/s 
which is consistent with the previous work^^. However as the sliding 
flake in our simulation is much larger than the nanoparticle (—250 
Au atoms) studied in the previous work^^, the Brownian motion can't 
activate the sliding flake overcoming the energy barrier of the com- 
mensurate state. As shown in Fig. 2, the kinetic energy in y direction 
and rotation reaches the same level when the sliding speed in x 
direction is smaller than 100 m/s. Therefore, the kinetic energy 
transfer among the three modes may cause the increasing of the 
sliding speed in x direction. 

Here, the intervals of superlubricity, Ts, vary widely. Detailed ana- 
lysis, summarized in Fig. 4b, shows that Ts depends primarily on the 
flake angular velocity, co, according to a simple relationship Ts = 60°/ 
cj. In other words, the rate of frictional scattering is simply deter- 
mined by the rate of crystallographic alignments. Roughly, Ts 
increases as the size of the sliding graphene flake increases. This is 
because the rotational inertia increases with the size of the graphene 
flake and for larger size flake the angular velocity is smaller. However 
the changes of the sliding speed are very complicated and highly 
dependent on the instantaneous state of the shding flake which will 
be systematically studied in our future works. 

The sensitivity of the results to initial temperature and flake ori- 
entation are given in Fig. 5. The simulated system is first equilibrated 
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The lower inset is configuration of nanoribbon substrate. The upper inset 
shows a small oscillation of the forward speed, (b) The rotation angle of the 
graphene flake for nanoribbon substrate and wide substrate. 

at the given temperature for 50ps by Nose-Hoover thermostat. Then 
the initial sliding speed of 400 m/s is added to the sliding flake and 
the thermostat is removed in the simulation afterward. The result for 
5 X Snm^andlO X 10 nm^ flakes launched initially at 0 K and 300 K 
show similar qualitative behavior, with regions of superlubricity 
punctuated by frictional scattering, see Fig. 5a. Detailed inspection 
shows a higher degree of noise for the high-temperature simulations, 
which is consistent with expectations for the impact of temperature. 
Simulations at T = 0 K, but for different initial orientations of the 
flake when it is launched, show overall similar qualitative behavior, 
with the sliding speed of all flakes decaying within 2ns, see Fig. 5b. 
For the initial misalignment between the sliding flake and the sub- 
strate (such as the rotation 30° in Fig. 5b) the initial superlubricity 
state can last a longer time. However after the first frictional scatter- 
ing it becomes the similar frictional scattering behavior as the initial 
alignment case. We have also tested other fiake shapes and varied 
substrate boundary conditions. In all cases we observe a similar 
stepwise decay, indicating that this is a computationally robust effect. 

Based on the above analysis, we conclude that frictional scattering 
is the main source of sliding friction of the graphene nanofiakes on 
graphite in the ballistic friction regime. Based on this insight, we now 
explore two routes to suppressing frictional scattering in order to 
achieve persistent superlubricity at high speed. 

Firstly, we introduce a graphene nanoribbon as a substrate to 
constrain the rotational motion of the graphene fiake during sliding. 



SCIENTIFIC REPORTS | 4 : 4875 | DOI: 1 0.1 038/srep04875 



5 




strain (%) 

Figure 7 | The stretching substrate, (a) The comparison of speed 
component versus time for stretched substrate and unstretched 
substrate. The lower inset is the configuration of 5% stretched substrate, 
(b) The shding potential for different stretched substrate, (c) The relation 
between the sliding energy barrier and stretched strain. 

Such graphene nanoribbons can be generated by the lithographic 
patterning^^ or chemical methods^^ with desired shapes and dimen- 
sions. In Fig. 6a, we launch a 10 X 10 nm^ flake on a graphitic 
nanoribbon of width 11 nm, with the flake rotated 90° relative to 
the normal AB stacking of the nanoribbon. The sliding speed of the 
flake decays by less than 3% with the temperature increasing less than 
1 K in the simulation time of 3 ns, compared with complete decay of 
the forward motion for the unconstrained flake. The nanoribbon acts 



as a frictional waveguide, periodically reflecting the flake off its edges, 
producing small velocity oscillations shown in the upper inset of 
Fig. 6a, and limiting flake rotation to ±1°, see Fig. 6b. However, 
for wide substrate as it can't constrain the rotation of the flake, the 
fictional scattering occurs when the flake first rotates 30° to encoun- 
ter the commensurate state, see Fig. 6b. Therefore the sliding velocity 
rapidly decreases as the results in Fig. 5b. The waveguide has a speed 
cutoff: smooth sliding ceases below —15 m/s, a typical velocity for 
transition to stick-slip motion in atomic-scale simulations^^. 

In a similar manner, we can achieve suppression of friction by 
biaxial stretching of the substrate. Such stretching of a graphitic 
surface can be achieved by uniformly stretching a PDMS film with 
graphene attached on it^\ For example, Kim et al.^^ have found the 
graphene can sustain as much as 30% uniaxial tensile strain without 
failure. Besides, as the lattice mismatch between the substrate and the 
graphene there is some residual strain in the epitaxial grown grahene 
which may be another routine to generate stretching in the graphene 
substrate^^. Because of the very weak interlayer van der Waals inter- 
action, the deformation of the substrate can not transfer to a flake 
adsorbed. So in this way we can envisage generating an initial lattice- 
constant mismatch between the substrate and the flake, as indicated 
by the moire patterns in the inset of Fig. 7a. We have modeled the 
sliding behavior of a 10 X 10 nm^ flake on a strained graphite sub- 
strate, with strains of 3% and 5%. For the 3% extension, the flake still 
undergoes scattering events, but with a reduced impact on the for- 
ward speed. For the 5% tension, there is no evidence of scattering 
events over the time of the simulation, 2ns. This is consistent with the 
notion that the flake can no longer be in exact commensuration with 
the substrate, due to its significantly different lattice parameter. The 
sliding potential £vdw for different stretched substrates is given in 
Fig. 7b. Here the orientation of the flake is same as the substrate and 
the sliding direction is in armchair direction same as the direction of 
the initial launched speed. By applying strain to the substrate the 
sliding potential barrier of the graphene flake is significantly reduced. 
As in Fig. 7c the sliding energy barrier for the 5% stretched substrate 
is only 5% of the value of unstretched substrate. This is the origin of 
low sliding resistance for stretched substrate. 

In conclusion, the motion of a nanoscale graphene flake on graph- 
ite at high speed is dominated by frictional scattering, where large 
oscillatory forces act each time the flake rotates through crystal- 
lographic alignment with the substrate. Frictional scattering is 
intrinsically two dimensional, and can be suppressed by constraining 
the flake to move along a quasi-one-dimensional frictional wave- 
guide or stretching the substrate to generate lattice mismatch. In 
addition, we have introduced a simple shearing mechanism as a 
practical way to accelerate nanoscale flakes of graphene to velocities 
> 100 m/s where ballistic friction occurs. These results should stimu- 
late further experimental and theoretical interest in achieving per- 
sistent superlubricity on the nanoscale, a regime of considerable 
practical importance for future nanomechanical devices. 

1. Urbakh, M., Klafter, J., Gourdon, D. & Israelachvili, J. The nonlinear nature of 
friction. Nature 430, 525-528 (2004). 

2. Luan, B. & Robbins, M. O. The breakdown of continuum models for mechanical 
contacts. Nature 435, 929-932 (2005). 

3. Mo, Y., Turner, K. T. & Szlufarska, I. Friction laws at the nanoscale. Nature 457, 
1116-1119 (2009). 

4. Socoliuc, A. et al. Atomic-scale control of friction by actuation of nanometer- sized 
contacts. Science 313, 207-210 (2006). 

5. Lantz, M. A., Wiesmann, D. & Gotsmann, B. Dynamic superlubricity and the 
elimination of wear on the nanoscale. Nature Nanotech. 4, 586-591 (2009). 

6. Hirano, M. [Superlubricity of Clean Surfaces] Superlubricity [A. Erdemir & J. M. 
Martin(eds.)] (Elsevier, Amsterdam, 2007). 

7. Shinjo, K. & Hirano, M. Dynamics of Friction - Superlubric State. Surf. Sci. 283, 
473-478 (1993). 

8. Hirano, M., Shinjo, K., Kaneko, R. & Murata, Y. Observation of superlubricity by 
scanning tunneling microscopy. Phys. Rev. Lett. 78, 1448-1451 (1997). 

9. Dienwiebel, M. et al. Superlubricity of graphite. Phys. Rev. Lett. 92, 126101 (2004). 



SCIENTIFIC REPORT: | 4 : 4875 | DOI: 1 0.1 038/srep04875 



6 



10. Depondt, P., Ghazali, A. & Levy, J. Self-locking of a modulated single overlayer in a 
nanotribology simulation. Surf. Sci. 419, 29-37 (1998). 

11. Filippov, A. E., Dienwiebel, M., Frenken, J. W. M., Klafter, J. & Urbakh, M. Torque 
and twist against superlubricity. Phys. Rev. Lett. 100, 046102 (2008). 

12. de Wijn, A. S., Fusco, C. & Fasolino, A. Stability of superlubric sliding on graphite. 
Phys. Rev. E 81, 046105 (2010). 

13. de Wijn, A. S., Fasolino, A., Filippov, A. & Urbakh, M. Low friction and rotational 
dynamics of crystalline flakes in solid lubrication. Europhys. Lett. 95, 66002 
(2011). 

14. Yang, J. et al. Observation of high-speed microscale superlubricity in graphite. 
Phys. Rev. Lett. 110, 255504 (2013). 

15. Liu, Z. et al. Observation of microscale superlubricity in graphite. Phys. Rev. Lett. 
108, 205503 (2012). 

16. Guerra, R., Tartaglino, U., Vanossi, A. & Tosatti, E. Ballistic nanofriction. Nature 
Mater. 9, 634-637 (2010). 

17. Lebedeva, 1. V. et al. Fast diffusion of a graphene flake on a graphene layer. Phys. 
Rev. B 82, 155460 (2010). 

18. Zheng, Q. S. et al. Self- retracting motion of graphite microflakes. Phys. Rev. Lett. 
100, 067205 (2008). 

19. Lindahl, E., Hess, B. & van der Spoel, D. GROMACS 3.0: a package for molecular 
simulation and trajectory analysis. /. Mol. Model. 7, 306-317 (2001). 

20. Mayo, S. L., Olafson, B. D. & Goddard, W. A. DREIDING: a generic force field for 
molecular simulations. /. Phys. Chem. 94, 8897-8909 (1990). 

21. Xu, Z. P. etal. Trans-phonon effects in ultra-fast nanodevices. Nanotechnology 19, 
255705 (2008). 

22. Guo, Y. J., Karasawa, N. & Goddard, W. A. Prediction of Fullerene Packing in C60 
and C70 Crystals. Nature 351, 464-467 (1991). 

23. Rivera, J. L., McCabe, C. & Cummings, P. T. Oscillatory behavior of double- walled 
nanotubes under extension: a simple nanoscale damped spring. Nano Lett. 3, 
1001-1005 (2003). 

24. Hedgeland, H. et al. Measurement of single- molecule frictional dissipation in a 
prototypical nanoscale system. Nature Phys. 5, 561-564 (2009). 

25. Woan, G. The Cambridge handbook of physics formulas. (Cambridge University 
Press, Cambridge, 2000). 

26. Miiser, M. H. Structural lubricity: Role of dimension and symmetry. Europhys. 
Lett. 66, 97 (2004). 

27. Miiser, M. H., Wenning, L. & Robbins, M. O. Simple microscopic theory of 
Amontons's laws for static friction. Phys. Rev. Lett. 86, 1295 (2001). 

28. Tapaszto, L., Dobrik, G., Lambin, P. & Biro, L. P. Tailoring the atomic structure of 
graphene nanoribbons by scanning tunnelling microscope lithography. Nature 
Nanotech. 3, 397-401 (2008). 



29. Li, X., Wang, X., Zhang, L., Lee, S. & Dai, H. Chemically derived, ultrasmooth 
graphene nanoribbon semiconductors. Science 319, 1229-1232 (2008). 

30. Braun, O., Peyrard, M., Bortolani, V., Franchini, A. & Vanossi, A. Transition from 
smooth sliding to stick- slip motion in a single frictional contact. Phys. Rev. E 72, 
056116 (2005). 

31. Kim, K. S. et al. Large-scale pattern growth of graphene films for stretchable 
transparent electrodes. Nature 457, 706-710 (2009). 

32. Pan, Y. et al. Highly Ordered, Millimeter- Scale, Continuous, Single- Crystalline 
Graphene Monolayer Formed on Ru (0001). Adv. Mater. 21, 2777-2780 (2009). 



Acknowledgments 

Y.L. acknowledges the support by the National Natural Science Foundation of China 
through Grants No. 11302163 and No. 11321062. Q.S.Z. acknowledges the financial 
support from National Natural Science Foundation of China through Grants No. 10832005, 
the National Basic Research Program of China through Grant No. 2013CB934200 and No. 
2007CB936803. 



Author contributions 

Y.L. conducted all of the MD simulations and prepared the figures. F.G. developed the 
frictional scattering model used to explain the results and wrote the main manuscript text. 
Q.S. supervised the study and wrote the manuscript. All authors reviewed the manuscript. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Liu, Y.L., Grey, F. & Zheng, Q.S. The high-speed sliding friction of 
graphene and novel routes to persistent superlubricity. Sci. Rep. 4, 4875; DDI: 10. 1038/ 
srep04875 (2014). 



©0®© 



This work is licensed under a Creative Commons Attribution-NonCommercial- 
NoDerivs 3.0 Unported License. The images in this article are included in the 
article's Creative Commons license, unless indicated otherwise in the image credit; 
if the image is not included under the Creative Commons license, users will need to 
obtain permission from the license holder in order to reproduce the image. To view 
a copy of this license, visit http://creativecommons.Org/licenses/by-nc-nd/3.0/ 



SCIENTIFIC REPORT; | 4 : 4875 | DOI: 1 0. 1 038/srep04875 



7 



